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AN EXPERIMENTAL INVESTIGATION 0? THE GAINS OF DIELECTRIC 

ANTENNAS 



INTRODUCTION 

The idea of using dielectric materials for the 
control of the direction of electromagnetic radiation 
was probably first conceived by Hertz in the Nine- 
teenth Century 1 , when he observed the strong influence 
of pitch blocks on this little understood form of 
energy flow. Subsequently, in the early part of the 
Twentieth Century a few investigations were reported 
concerning the use of dielectrics as wave guides, 
but there were only occasional contributions to this 
theory. Rapid progress In this direction was not 
begun until the availability of microwave equipment 
in the 1930's. By 1939, dielectric wave guide theory 
had been thoroughly studied and a few experiments 
had been attempted in the use of dielectrics as 
directional radiators. With the stimulus of the war, 
and the rapidly increasing need for more highly 
directive antennas, dielectric rods were studied, 



1. Superscript numbers refer to references listen in 
Bibliography. 
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tried, and put into use beginning in the early 1940’ s. 
In 1943 the U.S. Navy put into service a radar 
equipment, designed by the Bell Telephone Laboratories, 

l 

which employed an array of forty-two rods of 

Polystyrene to give a sharp directional beam in the 

horizontal plane. (Radar Equipment Mark VIII). 

The Germans, about the same time, began to use a 

2 

directi onal antenna on aircraft which consisted of 
four dielectric rods arranged to give a symmetrical 
pattern. Subsequently, a large number of dielectric 
antennas have been built, and an ever increasing 
demand for these small size, highly directional 
devices has encouraged more and more attention to 
their theory of operation. 

It should be pointed out in the beginning that 
^ the classical methods of antenna design, based on 
Maxwell's equations, have not been successfully 
applied to dielectric antennas whose dimensions have 
been of the order of a wavelength. As clearly pointed 
out by Stratton"^ in his book, "Electromagnetic Theory", 
Maxwell's equations are valid, in general, only to 
large scale phenomena, and attempts to extend them 
to the so called "Quasl-mlcroscoplc state" lead only 
to partial success. Although many of the laws of 
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classical electrodynamics apply to submicroscopic 
domains, it is not possible to apply them generally. 

On the other side, there is. quantum electrodynamics, 
which at this writing, have not been extended 
through the transition areas. The statistical 
averages, over a large number of atoms, of the 
equations of quantum electrodynamics ultimately 
must lead to the relations derived by Maxwell. 

Design of dielectric antennas to date has been 
based largely on cut and try methods. Some of their 
characteristics have been determined by various 
experimenters, and considerable theoretical work 
has been done in an endeavor to explain their 
operation. But there is no mass of data for the 
theorist to draw upon, nor is there an adeouate theory 
at the present time to enable calculation or prediction 
of the characteristics of a given dielectric shape. 

It is the purpose of the writer, in undertaking this 
work, to determine experimentally the ’’concentrating 
action” of certain of these antennas, and to express 
this in terms of gain relative to an isotropic 
oscillator; or, ’.. r hat amounts to the same thing, to 
express this gain as a ratio of effective area to 
actual area. It is believed that these data, thus 
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obtained, may be of use in subsequent investigations, 
or for other investigators, to ultimately develop 
a satisfactory theory of the function of the dielectric 
antenna. Finally, available theoretical approaches 
will be discussed objectively in the light of the 
experimental results obtained herein. 

V/ork of other investigators : 

zr 

Mallach , in Germany, presented a thesis in 
1939 which reported considerable work in the use of 
dielectrics to produce a directional beam. He also 
produced a simple theory which explains q uali tail velyy 
the action of the dielectric, but failed quantitatively 
to verify all of the known data. Taking a dielectric 
filled wave guide, v/hoee walls are considered 
perfectly conducting, excited at one end and open at 
the other, it is well known that radiation is emitted 
from the open end, and this radiation is directional 
along the projected axis of the guide. If the 
conducting sheath is removed from a part of the guide, 
nearest the open end, leaving some of the dielectric 
rod protruding, radiation will occur both from the 
end of the dielectric and from its side walls as 
well. Mallach^ theory is based on the consideration 
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that each element of the surface acts as a radiator 
and the phase of this excitation depends on the 
extended length of the rod.. He assumes in this work 
that the wave is linearly -solarized and of constant 
intensity and is transmitted down the rod without 
attenuation or reflection. Furthermore, he does not 
consider the rod cross section area except for the 
determi nation of the wavelength. He also considers 
only small angles of diffraction, taking no account 
of the spherical nature of the v/ave. Mallach 
reported a considerable amount of experimental data 
on both rods and dielectric tubes, and arrived at 
some simple formulas for choosing dimensions. For 
the sizes of rods and tubes that he used, his simple 
theory and formulas are satisfactory, but are of no 
use in obtaining lobe patterns for the entire 
radiation field nor for computing gains. 

• SouthworthT, at the Bell Telephone Laboratories, 
whose early publications on v/ave guides and v/ave 
guide theory are well known, has done considerable 
experimenting for more than a decade with tapered 
dielectric rods, both rectangular and circular, of 
diameters of a half wavelength or less. He and his 
associates have experimentally established that 
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the gain decreases as the attenuation in the rod 
increases. 

O 

Wilkes' , in a thesis presented to the University 
of Paris in 1946, examined the problem from the 
point of view of the lens action of the dielectric 
on radiation of wavelengths comparable to the 
dimensions of the lens. He pointed out that the 
apparent index of refraction plays an important 
role in the energy gain of the dielectric radiator. 

By plotting experimentally obtained energy gains 
versus length of exposed dielectric he has shown a 
correspondence of periodic gain to Snell's law of 
refraction. Furthermore he has shown that the 
apparent index of refraction is itself not constant, 
but depends on the exposed length of dielectric and 
a cross section dimension. A formula for calculating 
the apparent index from the real index of refraction 
is developed in his paper. Wilkes also computed some 
radiation patterns based on analogy to more 
conventional antennas, which verified very nicely the 
experimental data for the major lobe. But, as his 
paper points out, his assumptions used would cause 
considerable errors in the side lobes. This 



calculation assumes constant wavelength, which he 
shows to vary somewhat throughout the lens length. 

He also assumes, for this calculation, constant 
displacement currents. The complexity of an exact 
formulation of wavelength and intensity would lead 
to expressions which would be extremely difficult 

if not impossible to solve. 

IX X2 16 

Several investigators ’ ’ have subsequently 

published attempts at calculation of radiation 

patterns. All of the recent attempts known to the 

author are based on the formal solution of the field 

equations given in Schelkunof f * s well known book^, 

as applied to waves in dielectric wires, and expressed 

g 

in terms of cylinder functions. Elsasser in 1949 
published an analysis of the various non-radiatlve 
modes in a circular dielectric rod, in which he also 
pointed out a perturbation method of calculating 
dielectric loss. Elsasser, in this paper, did not 
attempt calculation of radiation from a finite length 
rod. This theoretical study is in substantial 
agreement with an experimental investigation of 
dielectric cylinders as wave guides performed by 
Chandler^ at RCA laboratories. 
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Marston^ of the U.S. Naval Research Laboratory 

has attempted calculation of radiation from dielectric 

rods of uniform circular cross section. He employed 

the method of assumed fields based on Schelkunoff ' s 

18 

Equivalence theorems , published in the Bell System 
Technical Journal in 1 936. Marston computed the 
radiation field arising from equivalent sources at 
the sides of the rod, and separately for those at 
the ends of the rod, and points out that maximum 
gain in the axial direction is achieved when the 
total radiation from the end is in phase with the 
radiation from sides (for modes of type n = 1 only). 
Marston’ s calculations also show that the actual 
values of gain for dielectric circular rods cannot 
in general be predicted from ordinary array theory, 
but prospects of achieving this exist through the use 
of his complex expressions for the components of 
electric field intensity. However, all of Marston* s 
assumptions are not in agreement with findings of 
other investigators; for example he assumes constant 
wavelength within the dielectric body. If this 
variable were substituted for the constant, it would 
make his equations almost hopelessly complex. 
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A great deal of the known Information on 

12 

dielectric antennas has been assembled by McKinney 
in a dissertation presented to the University of 
Texas in June, 1950. He reported also a very 
substantial a.mount of experimental data, as well as 
a slightly different method of attempted calculation 
of the fields. The method presented is also based 
on Schelkunof f ’ s Equivalence Theorems, but uses an 
extension devised by Horton, in which he computes 
fictitious electric and magnetic current sheets 
about the dielectric which are a useful aid in the 
analysis. Then he obtains the vector potentials at 
a point in space due to these sheets and computes the 
E and H fields at the point. In his calculations he 
neglects the radiation from the end of the rod, 
neglects reflections in the rod, and neglects any 
changes which might occur to the wavelength of the 
wave within the dielectric. McKinney's experimental 
data agreed with his theoretical calculations quite 
closely. For the sizes and shapes used and for the 
TMqi mode, it appears that his solutions are quite 
precise. However, his report does not show any 
calculations of patterns for dielectric antennas 
excited in the dominant mode, which is the mode most 
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commonly used, nor does his theory contain a precise 
statement of the optimum gains which may be obtained. 
The TMq^ mode gives a null on axis, which is 
unacceptable for most cases where dielectric antennas 
are required. 

It should be pointed out that no theoretical 
treatment of any dielectric antennas other than those 
of circular cross section have been found in any of 
the literature. As will be brought out later in this 
essay, an attempt at a straightf orward solution to 
the rectangular shaped antenna by any of the classical 
methods would lead to exceedingly complex equations 
whose solutions would involve serious difficulties. 
Nevertheless, in a great many applications, the 
rectangular antenna is the one preferred due to the 
predominant use of rectangular wave guide which 
assures no rotation of the plane of polarization. 

Hence this essay includes experimental gains obtained 
for both rectangular and circular cross section 
guides, including such correlation to existing 
theories as are possible. 
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CHAPTER 1. 



REVIEW OF BASIC LAV’S OF ELECTROMAGNETIC RADIATION 



Th© objective of Chapter 1 Is to develop the 
equations for the propagation of electromagnetic 
waves in a circular dielectric rod, infinitely 
long, imbedded in a homogeneous medium. The starting 
point will be Maxwell's equations. The development 
will be carried out for the cylindrical case because 
treatments of this are available in the literature. 
Throughout this work the MKS system of units will 
be used. 

It has been pointed out in the introduction 
that Maxwell's equations hold with certainty only 
for macroscopic phenomena. Accordingly, we make 

the assumption throughout Chapters 1 and 2 that only 

/ 

macroscopic conditions are considered, that the 
media considered are isotropic and homogeneous, 
and that, except for stated discontinuities, the 
physical properties in the neighborhood of every 
point are continuous. Under these conditions, the 
following Maxwell's equations have been firmly 
established. 



( 1 ) 



V*E 




t 



( 2 ) 



V* H - = J 
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whore the vector quantities E, B, H, D and J are 
respectively the electric intensity, the magnetic 
induction, the magnetic intensity, the electric 
displacement and the current density. 

Since the divergence of the curl of any vector 
is identically zero, we obtain by taking the 
divergence of (l) 



(3) v- = o 



and under the condition of continuity imposed above 

(4) 



dt 



V • B - 0 



Since the divergence of the magnetio intensity does 
not change with time, and certainly at some time 
it did not exist, it follows that 

(5) v- 8 =o 

Likewise, taking the divergence of (2) and 
commuting the partial and del signs 

l6> V- J- + -t? V- D = ° 



15 



Nov/ If we denote the charge density at any 
point by f, the equation of continuity or the 
conservation of charge states that the time rate of 
change of Pplus the divergence of the current density 
is identically zero, or 

(7) V • xJ + = 0 

Substituting for V--J in (6) 

h v-d - U = o 

or 

(9) f-j (V-D'fj = 0 

Applying again the argument that if the time rate 
of change of the quantity is zero, the quantity 
itself must be zero if it ever has vanished, it 
follows that 

do) v- D - ? 

This tells us that the charges distributed with the 
density f constitute the sources of the vector D. 

51 ec trio and Magnetic Polarl zatlon 

In order to describe the electromagnetic state 
of a substance, it is convenient to define two 
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additional vectors, the electric and the magnetic 
polarisation vectors. By defini ion 

(11) P - 0 - £c E 

da) m - 4 e -h 

where £ a is the inductive capacity in free space 
and M 0 is the permeability of free space. Therefore 
in free space, the polarization vectors vanish. 

Substituting the values for D and H from (ll) 
and (12) into (2) and (10), vie obtain the following 

( 13 ) V* §£■ = m. (j > §% +Vx/m) 

(14) v • e - - v- p; 

In any physical media which is isotropic and 
in which the field vectors are linear, the Electric 
Displacement and the Magnetic Induction are related 
to the corresponding field intensity by the constants 
of the medium 



( 15 ) 



D = £E 
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( 16 ) B - m H 

From the above equations (13-14) it is possible 
to account for material bodies in an electromagnetic 
field by -an equivalent charge of density -VP and 
an equivalent distribution of current density -§£ t 
It is convenient at this point to define two 
additional coefficients of proportionality. The 
electric and magnetic susceptibilities of any media, 
excluding ferromagnetic material, will be denoted 
by x Q andl^ and are defined by 

dr) ?-- x e e.E 

(18) M 

We can define these scalar quantities in this way, 
since in isotropic media, the polarization vectors 
are parallel and proportional to the corresponding 
fifeld vectors. 

It is also convenient to define the relative 
dielectric constant 

(19) 6 = £ 

and the relative permeability 



( 20 ) 



V - 
"to ~ Mo 
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From (ll), (15), (17), and (19) It follows that 

( 21 ) X e = £ -I 

and from (12), (16), (18), and (20) that 

( 28 ) X„ - K„- I 



Vector and Scalar Potentials 

To facilitate the analysis of electromagnetic 
fields it is useful to define some auxiliary functions 
called potentials. These may be shown as follows: 
Equation (5) shows that the divergence of the 
Magnetic Induction Vector, B, is zero at all times. 
According to this, the field of B is always 3olenoidal, 
and therefore B can be represented as the curl of 
another vector A. 

(23) B = VXA 

Substituting this value of B into equation (l) 

gives 

(24) V* (E+ §jf) - 0 

Thus the field of the vector (E f -|j ) is 



irrotational and equal to the gradient of a 



scalar 
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function which we may call <f). 

(25) 2 f = - V ^ or E = -V^-§£ 

The function A is called the vector potential of 
the field, and the (f) is the scalar potential. 

It should be noted, however, that the vector A 
and the scalar (f> are not unique. Thus B could have 
been represented by the curl of another vector 

(26) B = V X A 0 

where 

(27) k=A 0 -Vf 

where ^ is any arbitrary scalar function of position, 
It would similarly follow that 

(28) E * - V <f> 0 - ^ 

and <f) and (j) Q would be related by 

(29) <t> = f If. 

Likewise an infinite number of potentials describing 
the same field could be derived in the same manner. 

By the definition of<f, equations (15), and 
using the above defined potential, equation (25), 
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(jo) s> = - e(vtj> + If) 

And by equations (16) and ( 23 ) 

(31) H - yc V XA 

Substituting (30) and ( 31 ) into eouation ( 2 ) 
and multiplying the above by/4 and commuting 

Cr 

v/ith V gives 

( 32 ) V*W + -«£vf| M J 

By substituting ( 30 ) into equation (10) and dividing 
by - 6 

(33) V- \7i> f V-ji - ? 



Nov; all particular solutions of ( 32 ) and ( 33 ) 
lead to the same electromagnetic field v;hen subjected 
to the same boundary conditions. They differ among 
themselves only by the arbitrary function . In order 
to specify the functions more precisely, the 
condition may be imposed that A and also satisfy 
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( 34 ) v-A t^E =0 

In order to do this, it is necessary that <p satisfy 

(35) vy -MSrjr> = V-Ac 

Where A Q and <fi 0 are particular solutions of (32) 
and (33). 

The potentials A and (f> hy the aid of the above 
are completely specified. They are solutions of 
the following equations, which are the same as (32) 
and (33) except for the substitution of 
from (34) . 

(36) V*V*/I - VV-A+*£$j: 

( 37 ) vy -me 0 =-i? 



3y use of the vector identity 

(38) V*VXA = VV-A -V-VA 

equation ( 36 ) reduces to 

1 

(39) V'A-Mt^j* = -M-J 
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The relations (23) and (25) for the vectors 
B and E are not general. To them may be added any 
particular solution of the homogeneous equations: 

vxe : t §7 = 0 

(41) Vx H - - 0 

(42) V - B = 0 

(*3) V- p = 0 

From the symmetry of this system it is evident that 
it can be satisfied identically by 

(44) d = - S 7 M* 

(45) H = - 

>■) 

from which may be obtained 

(46) E = - -g- 



(47) 



B = 
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But the new potentials, A* and ( p must be subject to 
the conditions 

(48) -JC& fjp- ~ 0 

(49) -Mi jrp. = 0 

(50) V- A* 

A general solution of the inhomogeneous system 
given by equations (l), (2), (5) and (10) is, then 

(51) B* VAA - M -§y - M V<f>* 

(52) t E=-V^- §f - ■£ VM* 

The functions A - * and (f) * are potentials of a source 
distribution entirely external to the region 
considered. Usually A and (p are set equal to 
zero, and the potentials of all charges, wherever the 
source, are represented by A and <j). 

At any ooint where the charge and current 
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densities are zero, a possible combination of 
is , 

potentials^ A Q = 0 and (p 0 = 0.* The function^ is 

then according to (35) any solution of the homogeneous 

equation 

( 53 ) vy - Mgffr = o 

Since (f> Q = 0 satisfies this same equation at this 
pointy may be chosen in such a way that <f) vs,nishes. 

In this case the field can be expressed in terms of 
only the vector potential. Equations (23) and (25) 
become 



(54) 


b = V*A 


(55) 


„ 14 

2 = “ 3+ 



and (39) and (34) become 

(56) 

(57) V-A - 0 

The foregoing paragraphs have shown how an electro- 
magnetic field' in a homogeneous medium may be 
reduced to a vector and a scalar ootential which 
will satisfy Maxwell’s eauations. 
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Hertz Vectors or Polarl zatton Potential s 

It is possible under ordinary conditions to 
define an electromagnetic field in a homogeneous 
medium by means of a vector function alone. 

Considering at first, no free charges or 
conduction currents v/ithin the medium, the field 
equations reduce to the homogeneous form (equations 
(4o) through (43) ). Assume a vector, II, whose 
time derivative is proportional to the vector 
potential A 

(58) A = m£ §#• 
and substitute in (23) and (25) 

(59) B * Mi V/ ff 

(60) E= -V$-m£xP- 



Using the relations B = -tt//and D = iE, substitute the 
above values in the homogeneous equation (4l). After 
commuting the and V signs, the result is 

(61) If (V* V*U t Vtj) = 0 
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Since the scalar function <f) is arbitrary as long 
as it satisfies equation ( 34 ), it may now be chosen 
to also satisfy 

( 62 ) (j, = - V-Z 

If the above chosen value of (ft is substituted in 
equation (61) and the result integrated with respect 
to time, the following is obtained 

'^jr 

(63) vx. VAJ - 7V'I ^Ti - a constant 

A 

The value of the constant does not affect the 
determination of the field, and so it may be replaced 
by zero. Every solution of the vector equation 

(64) V*V*jr - V V‘E 0*-- 0 

determined an electric field through 

(65) B = V*# 



and 

(66) E= | ^ 



Equation ( 65 ) Is identical to (59) and (66) is the 
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same as (60) with the substitution of VV'JT for 
its equal -70 by virtue of (62). 

Applying the vector identity, equation (38), to 
equation (64) gives 

(67) V- VI -MS ffi = 0 

The requirement that ^ satisfy equation (53) is 
met in equation (64) and (67). 

Inasmuch as the Electric Displacement Vector, 
D, as well as the Magnetic Induction is solenoidal 
in a charge free region, an* alternative solution 
could be similarly constructed of the form 

(68) A” - mB | vf 

(69) f « - V-r* 

(70) 

(71) H« VV 

where II 



is any solution of (64) or ( 67 ). 
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It therefore follows that the electromagnetic 
field throughout a homogeneous Isotropic region in 
which no charge nor current exist could he resolved 
into two partial fields, one derived from the 
vector II and the other from II ' r . The source of 
these fields is exterior to the region. In order to 
show the physical significance of the Hertz vectors 
it is necessary to relate them to their sources. 

This may be done by finding inhomogeneous ecuations 
from which ( 67 ) is derived. 

The vector D may be expressed in terms of the 
electric field intensity and the electric polarization 
vectors, using equation (11), as D =<f 0 E + P. Then 
substitute this in the homogeneous field equation. 
Equation (41) becomes 

(72) v*/W„i£=££ 
and (43) becomes 

(73) V-E = -£ v-p 

The above two equations as well as (40) and (42) 
are still identically satisfied by (65) and (66), 
provided B Is replaced by £ and that II is now any 



0 
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solution of 

(74) 7' VI -*£.& =-*/> 

Tho source of the vector II and the electromagnetic 
field which it prescribes is a distribution of 
electric polarization P. 

The vector P is the electric dipole moment per unit 
volume of the medium. The vector II is associated 
with a distribution of electric dipoles and may be 
called the electric polarization potential. 

Similarly the field associated with II is 
established by a distribution of magnetic polarization. 
The magnetic induction vector, B, is related to the 
magnetic field intensity, H, according to equation? 
(12), by B - (H f K) . Substituting this Into 
(40) and (42) gives 

(75) 

and 

(76) V-H = - V-/M 

If it is further stipulated that II* shall be a 



solution of 
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(77) V' V I* T7* = ~M 

and that^ replaced in eouations (70) and ( 71 ) » 
o 

then equations (75), (76), (4l) and (43) are 
Identically satisfied by (70) and (71). 

The electric polarisation P may be induced in 
the dielectric by the field E, but it may also 
contain a part whose magnitude is governed by 
entirely external conditions. Insofar as dielectric 
antennas are concerned, only the Independent part, P Q , 
is of interest, which represents the electric moment 
of dipole oscillation activated by external power 
sources. This is also true of magnetic polarisation. 
To represent these conditions, equations (ll) and (12), 
expressed in modified form are 

(78) E=-^(D-r 0 ) 

(79) H = A B ‘Mo 

in which P and are ore scribed and independent 
o 0 

of E and H and where the induced polarizations of the 
medium have been absorbed into the parameters £ and 
M, . The electromagnetic field due to these distri- 
butions of P 0 and M 0 are given by 
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(80) 



E - 



vv-e-m , f 0 



~M V* 



d-f 



( Q1 ^ 



Equation ( 8 o) above was formed, by substituting in 
eouatlon (78) the value of D given by (70) and P Q 
given by (74), replacing,/^ and £ 0 by yU. and 
Equation (81) was formed by substituting into (79) 
the value of B given by ( 65 ) and M Q from ( 77 ). 

For these eauations to be valid, however, it must 
also be stipulated that II and II be solutions of 
equations corresponding to (74) 

(82) V' VI -y»e 0 =-£ P. 

(83) V' V1*-m£ §£ *-Mo 

Using the vector identity, equation ( 58 ), and 
substituting ( 83 ), it ±&salso possible to write ( 81 ) as 
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(84) H = £V + V'Vtll* -Alt, 



The last relation shows that, since B = VXA y the 
vector potential A can be derived from the Hertzian 
vectors by setting 



Boundary Conditions Discontinuities in the Field Vectors 

The above equations have all been based upon 
the assumption of continuous and isotropic media. 

It is now necessary to consider the relations which 
arise due to the boundaries between two media, where 
there occur abrupt differences in the inductive 
capacities, X and g , and in the conductivity, 6 . 

To aid in the analysis, it is useful to postulate 
a very thin transition layer, S, between the media. 



(85) ^ = Vxl* - VI// 



The associated scalar potential function is 



(86) <j> = - VJZ + |f 



in which ^ must also satisfy 



(87) 
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wherein the parameters^ , £ , and <s , vary rapidly but 
continuously from medium 1 to medium 2. Within this 
layer, S, the field vectors and their first derivatives 
are continuous, bounded functions of time and space. 
Consider a small right cylinder, Figure 1, in the 
layer such that its end surfaces lie in the surfaces 
of the layer and separated by a length A £ . 

From the volume integral of eouation (5) and 

IQ 

the divergence theorem of vector analysis , it 
follows that 



where the second integral is a surface integral 
over the walls and the ends of the cylinder, H being 
the unit vector normal to S' . If the elemental 
area of the end of the cylinder is sufficiently small, 
B may be considered to have a constant value over 
the end. Neglecting differentials of higher order, 

( 88) may be approximated by 



( 88 ) 




(89) ( B * n, + 6 Ht) A C{ + contribution of walls = 0 



F/q / 




Normal Boundary Cond/t/otf 
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where f7 ( Is unit vector normal to the surface In 
the sense inward to medium 1, and \) z the unit 
vector normal in the sense inward to medium 2. 

Since the contribution of the walls is proportional 
to , it will become vanishingly small when the 

limit is taken -» 0 , and the transition layer 
shrinks to the surface S. Designating the value 
of magnetic induction at the point on S in medium 1 
as B^, and in medium 2 as B^, and passing to the 
limit 0 and A&-+Q 

(90) (B a - BJ • n = 0 

Thus the transition of the normal component of B 
across any surface of discontinuity in the medium 
is continuous. 

Applying the same process of reasoning to the 
electric displacement vector D, equation ( 10) is 
integrated over a volume and its equivalent from 
the divergence theorem substituted. In this case, 
however, since V- D is equal to the charge density^ P , 
when it is integrated over 'a volume it becomes 
equal to the total charge within the volume, q, thus 

(91) f^VDd/r - Ion da 
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This charge, a, is distributed throughout the 
transition layer with a density p , and as the 
cylinder is reduced toward the limit, the q must 
remain the same because of the conservation of 
charge. Thus 

(92) ^ = f A 9 A a 

As a % 0 , f must then-**®. It is then convenient 

to replace the product, f aX , by a surface density 
or charge per unit area,^ . The transition of the 
normal component of the electric displacement vector, 
D, across any surface, S, is then given by 

(93) (D z - D.) n = w 

The presence of a layer of charge on S results in 
an abrupt change in the normal component of D, the 
amount of the discontinuity being equal to the 
surface density in coulombs per square meter. 

It is also necessary to investigate the 
behavior of the tangential components of E and H at 
the boundary. To aid this analysis, consider a 
rectangle of length AS and thickness aX > instead 
of a cylinder, lying at the surfaces of the two 
adjoining media. Let 5 0 be the a.rea of this 
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rectangle and C 0 be its closed contour. 

Taking equation (l) and integrating over the 
surface 



(94) f n 0 dc i +/ if-nda = 0 

Jo 'So 

where r? 0 is the unit normal vector in the conventional 

positive sense, as shown in Figure 2. By Stokes 

19 

theorem , this first surface integral is equal to 
a line integral of the vector E about a closed 
contour enclosing the surface. Thus (94) may be 
expressed as 

(95) f Eds + / f 7 -rt. da = 0 

Jc 0 /So 

Designate by y and 71 the unit vectors in the 
direction of circulation as shown in Figure 2. Then, 
neglecting higher order differentials, equation (95) 
may be expressed by 



(96) (£ y t contribution from ends = - §j aSaQ. 



Fig. B 



(?) Si Ut Si 



(!) S, M, 6, 




Tantgeniial Boundary C ond/fm 
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But the contributions from the ends are proportional 
to aX which becomes vanishingly small when the limit 
is taken. Let n be the unit vector normal to S 
as shown, and T be the unit tangent vector such that 

(97) r = n a *n 

Also, it is a vector Identity that 

(98) n o x n • E - n 0 • n * E 

Thus, when the limit is taken a$-> 0 , AS -*■ 0 

(99) n. ■[«'(£*-£,) + -0 



Since fl 0 ^ 0 and the orientation of the rectangle 
is arbitrary, the bracketed expression above must 
be zero, and so 

(100) n * (E 2 -E) = -Urn If- aH 

And since we have imposed, the conditions that the 
field vectors and their derivatives be bounded, the 
right side of the above equation vanishes when aX-*0 
Therefore 

(101) nx(E z -£,)-0 
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and it is seen that the tangential components of 
the vector S through a surface of discontinuity 
are continuous. 

The behavior of the magnetic intensity vector, 

H, at the boundary can be deduced by the same 
process. Integrating equation (2) over the surface 
of the rectangle described above and applying Stokes 
theorem to the first term, there is obtained 

( 102) j j-j- - J - J -J* H„ dtf 

- 5 © Jo 

and It similarly follows that 

( 103) n * (Wz -H.) = bm (§£ + j) aI 

The limit of -§7 A $ as Af~> 0 will vanish 
because D and its derivatives are bounded. But it is 
oossible that jj)T) -J a X may not vanish, because 

the current I = J-n 0 A 5 A$. is squeezed into an 
infinitesimal layer on the surface as the sides 
are brought together, and -J -» °° . It is convenient 
to represent this surface current density by K, 
which is defined as the limit of -Ja$ as 
approaches zero and J becomes infinite. Then 
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(io4) n*(H z -H) =K 



When the conductivities of the adjoining media 
are finite, no surface current can exist, because 
S is bounded and the product of 6EaI vanishes 
with aH . In the analysis of dielectric antennas, 
this condition exists, and so (104) can be reduced 
to 



In the foregoing, expressions were developed 
for the transition of the normal components of B 
and D, equations (90) and (93)» and for the 
tangential components of E and (H), equations (101) 
and (105). By using ( 15) and (16), the normal 
components of H and E are immediately obtained for 
the transition of an electromagnetic field from 
one medium to another. 



( 105) n K (H z ~ Hj -0 



( 106 ) 




(107) 
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and for the tangential components of D and B 



( 108 ) n x 



0>* - % P,) -0 



(109) 




e,) =o 



( for finite £ ) 



Poynting 1 s Theorem in regions free of Hysteresis 

Poyntlng’s theorem- gives the relation between 
the rate of change of energy stored in a field and 
the energy flow. This may be obtained as a general 
integral of the field equations (l), (2), (5), and 
( 10) . Since £ has the dimensions of mover 
expended per unit volume, it is apparent that this 
may be obtained by a scalar multiplication of (2) by 
S. Thus 



Likewise actuation (l) can be given the units of 
power per unit volume by multiplying it by H, Thus 



(no) E-V*H - E-|f = E-J 



(in) H- V * E + H • |f = 0 
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Subtracting the above two equations 

(112) H-V*£ -E- VaH +£-i -*■ H- If =-E-J 

But, by vector identity, the first two terms of 
(112) are equal to V-(E*h) and so this equation may 
be written as 

(113) V-(E*H) + E-J = - E-ff - H- 0 - 

Integrating over a volume V bounded by a surface 
S, and applying the divergence theorem to the first 
term, this becomes 



(ii4) fjE'H )-nda tfe- Jd* = -[(E-§ 

The terms on the right above represent the rate of 
decrease of the electric and magnetic energy stored 
in the volume. The second term on the left expresses 
the sum of the losses in Joule heat and the power 
expended by the flow of charge against the impressed 
forces. If there are no transformations of 
electromagnetic energy into elastic energy of a 
stressed medium, the balance is maintained by a flow 
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of electromagnetic energy across the surface bounding 
the volume. This is Poynting's interpretation of 
the significance of the surface integral in (114). 

The reduction of electromagnetic energy stored in 
the volume is partly accounted for by the Joule 
heat loss, partly by the energy introduced by the 
impressed forces, and the remainder flows outward 
across the bounding surface S, which is represented by 

( 115 ) j£ 5 -nda = l (E*H )-nd* 

Now the Poynting Vector is defined by 

(116) 6 r E x H watts/meter 2 

which may be interpreted as the Intensity of energy 
flow at a point in the field, or in other words, 
the energy per unit time crossing a unit area whose 
normal is oriented in the direction of the vector 
E X H . It is Poynting's theorem that the vector 
product (116) at any point is a measure of the rate 
of energy flow per unit area at that point. 

Another useful form of Poynting's theorem is 
the differential form, in which the Joule heat 
dissipation is separated from the power expended 
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by the flow of charge against the impressed forces. 
While E represents the total electric field 
intensity, let E be the Intensity of the impressed 
electromotive forces. Then 

(117) E * T -El' 

and 

(118) e-j=£-e'j 

Substituting (118) into (113) and replacing D by 
its equal ££and B by^/y, the differential form 
becomes 

(119) 

V'(E*h) + * f?(r H’) = EW 



The Complex Poyntlnx Vector 

Let h -M+w represent the density of the 
electromagnetic energy at any instant and ft -=g 
be the pov. r er expended per unit volume in Joule heat 
and by the flow of charge against impressed forces. 
From (119) with these substitutions, the Poynting 
theorem from ( 113 ) can be written 
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(120) V'5 -I- f£ t(p = o 



In a stationary field, h Is independent of time and 
the middle term would be absent. Moreover, when 
the field is a simple harmonic variation of time, 
the time average of the energy density over a period 



indicating an average value for a period. In the 
case of periodic fields, It follows that 



Now these field vectors are harmonic functions of 
time and can be represented In complex form with 
real and imaginary parts. By complex variable 
theory, it can be shown that the time average of 
the product of the real parts of two complex 
quantities is equal to one half the real part of 
the product of the first quantity times the 
conjugate of the second quantity. Symbolically 




V-5 + (? ^0 



(121) or 




( 122) fa (/I,) ■ fa (A z ) - i fa (A, A z ) 
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where is the comolex conjugate of Ag. 

According to the above, the mean intensi ty of 
the energy flow in a harmonic electromagnetic field 
is the real mart of the complex vector EXH 

<5 = (fto. (E)x 61a (tj) 

W23) „ 

S = jr a*(€* H) 

This is the complex Poynting vector which may also 
be denoted by 

(124) = 



In a plane wave the magnitude of H i 3 ~~ , where 
^ is the impedance of the medium. Therefore in 
engineering practice the expression 

\Z 



(125) 



Power 






is frequently used, corresponding to the vector 
expression (124). 
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Representation of a cylindrical field by Hertz Vectors 

In a homogeneous medium which Is free from 
any fixed polarization, the Maxwell Equations (1) 
and (2) are satisfied by the following relations In 
terras of Hertz vectors 

(126) £ = Vx V*! -A Vx If* 

(12?) H = Vx(£§f +el) 1 ■ V* VxiT*' 

The first of these can be readily verified by taking 
equation (80), applying the vector identity (38) 
to the first term, and noting that ~m£ 

by ( 67 ). The second one is similarly verified. 

Now consider a set of coordinate surfaces. 

Figure 3 » formed by a family of cylinders whose 

elements are parallel to the z axis. These 

cylindrical surfaces need not necessarily be 

circular or closed. Unit vectors i, , i„ , and i_ 

id 3 

are assigned as shown, i^ being normal to the 

z 

cylinder, i^ tangent to it and parallel to the/iaxis, 
and lg is tangent to the surface and perpendicular 



FIGURE 3. 




Relation of Unit Vectors to a Cylindrical Surface. 
The generators of tho surface are parallel to the 



£ axis 
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to l n and 1^ , forming a right handed orthogonal 

set. Position with respect to axes determined by 

the three vectors is measured by the coordinates 
1 o 

A , M~ , j? . An infinitesimal line element is 

(128) ds - U, M‘ + a«* ! 1 i, At 

Consider an electromagnetic field in an isotropic, 

homogeneous medium, unbounded in extent, existing 

such that the Hertz vector, II, is directed along 

the 2 axis. I^ = 0 II. =• 0, and II* - ^ . 

* 

II =0. The electric field intensity vector of 
this field, as given by (126) is 

(129) e ( “ = vxvxjr 

and the magnetic vector of the field, by (127) is 

(uo) H 0> = (£& *e) VxH 

The components of the vectors E can be calculated 
from (129) by means of the vector formula for the 
curl in an orthogonal system of curvilinear 
coordinates (reference 19, page 90 ), applying this 

twice. 
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( 131 ) 



vx f — \ t (fulfil — 2(^3 fc) \ i 63 f ^(SzJzL 

"&/?3V ^ ^3 A y^/ / 

^1 c ^ "&a = ^T,. ^ 




Since II, = U z ~ 0 , this comes out to be, after 

taking the curl twice 



cO) _ ± Wt 
1 JudiZu} 



( 132 ) 




I aVT z 




* 



Likev. r ise, to compute the components of H from (130), 
simply take the curl once and multiply by j 



thus 
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( 133 ) 

= 0 

From the z scalar comnonent of a function 

U z ~ */" » an electromagnetic field has been developed 

which is characterized by the absence of any axial 

or longitudinal component of the magnetic vector. 

It was earlier pointed out that II is the electric 

polarization potential, and therefore this field 

may appropriately be called a field of the electric 

4 

type, but following the terminology of Schellcunoff , 
it is more often referred to as the "transverse 
magnetic field." 

It is necessary to digress to develop here 
the wave equation, which the function K 2 = 
must satisfy. 

Taking the curl of equation (l) and substituting 
aH for B 



(134) Vt(VXE) = 
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Sine© J = &E and D -£E , equation ( 2 ) can be 

written 

U35) Vx H v£-§£ 

Substitute (135) into (134) and making use of the 
vector identity (40) 

( 136 ) v(v-e) - V- = -^<£T-§f 

Rewriting equation (10), substituting EE for D 

« 

(137) V •£ = 

Since in free space or in any conducting medium f 
is independent of the field distribution, the first 
term of ( 136 ) may be taken to be equal to zero. 

Then ( 136 ) becomes 

( 138 ) <7% -MS 4§r 

Since -J ~ sE » this could also be written 



(139) =0 



? 



4-9. 



By eliminating E from Maxwell's equations, there is 
also obtained 



(140) 




~0 



Now, returning to the equations of the cylindrical 
field, the components of E and H given in 
(132) and (133) must satisfy the scalar wave enuation 



(141) V> % -0 • 



The Laplacian of the scalar^, expressed in general 
orthogonal coordinates (reference 19, page ISO ),is 

(142) 




Substituting this in (l4l), where = 1, 0 ?- 2 

(143) 



_L ± 



Qu. ) 



+ 



> & ft 1 , 

Mi {%> *<y 



5? = 



0 



The elementary harmonic solutions of (143) 



are of the form 
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(144) ), = 



where {(u'jU*) 



is a solution of 



(145) 



1 

JtffodU 1 U, **y 



I J 3 f'fij \ , 



0 



Similarly, a partial field can be derived from a 

* 

second Hertz vector II by the operations 
(146) E® - -A(h V*l* 

H (a) = V A V X J 7 * 

If the vector II W is directed along the 2 axis, 
the vector components come out to be 



(147) 



_ _ 2L d*JT/ 

r 2£ ^ 

~ ^ a/eu' 




0 



( 148 ) 




i fi? 




i r i_ f Si 

J j x Ldu’[ a { 



dZ?) , d (*, Ml ) n 

d<SJ r W\ 4 z au*/ J 



The scalar function II* Is a solution of (l4l), 
and is a field of the magnetic type or transverse 
electric, characterized by a zero longitudinal 
component of E. 

If the partial fields derived from II ? and 11^ 
are superposed, the electromagnetic fields are of 
such generality that they can satisfy any prescribed 
set of boundary conditions on any cylindrical surface 
whose generating elements are parallel to the Z 
axis. In general, the choice of the coordinates is 
limited practically to systems in which the va-riables 
of equation (145) are separable. It is for this 
reason that this development of equations of electro- 
magnetic waves in a rod is limited to the case of 
the circular cylinder. 
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Wave Functions of the Circular Cylinder 

The simplest of the separable cases is that in 
which the family u' = constant is represented by a set 
of coaxial cylinders, where t< = r > = I 

and , See Figure 4. In this case equation 

(145) reduces to: 

(149) 

I f?(> ff) + -f fir* + (4 a -4f = 0 

which can be separated by writing f (V, &) as a product 

(150) f -fete) fete) 

where f t (r) is an arbitrary solution of the eouation 

(151) 

r w (r ifr) = 0 

and fe (&) is an arbitrary solution of 

(152) + Vfe = 0 

Solutions of (151) are Bessel’s functions, and 
a solution to (152) is fete) - 6 x ' Y ' e . The parameter, n, 
is a separation constant. Its choice is governed by 
the physical requirement that the field must be 
single valued at any fixed point in space. Since at 
the present, the discussion is limited to continuous 



FIGURE 4. 




Coordinates of the Circular Cylinder 
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media, the field is periodic in 0 and the values of 

n are limited to integers, 0, 1, 2, 3 either 

positive or negative. However, in a field represented 
by particular solutionsof (151) and (152), which 
is bounded by planes 6? -<9, , and Q - , non-integral 

values would in general have to be assigned to n. 

The constant h is, in general, a complex quantity 
pertaining to the propagation of the wave and is 
called the propagation factor. An explicit expression 
for h in terms of the frequency and the constants 
of the medium can be obtained only after '"he behavior 
of if/ over a cylinder constant on a plane ?= 
constant has been prescribed. 

The equation (151) is Bessel's equation, the 
solutions of which are the Bessel functions r) • 

In accordance with terminology enrol oyed by Jahnke 
and Emde and several other authors the term "circular 
cylinder function" will be used to denote any 
particular solution of (151) and designated by the 
letters fj - (^-J\ f) . The order of the functions 

whose argument is r is n. Particular solutions 

of the scalar wave equation (141) which are periodic 
in t and 0 can be constructed from elementary waves 



of the form 
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(153) 

e *”* 2*, (I CFR r) e tA * z 



Reference i 3 made to standard texts In 
mathematics for the properties of the various Bessel 
functions, but it is essential to this development 
to note several pertinent characteristics. The 
Bessel function, or circular cylinder function of 
the first kind, is finite and converges for all 
finite values of the argument. For non-integral 
values of n it has tvro linearly independent 
solutions for -^(f) and ; but if n is an 

Integer the tv;o solutions are not independent. 

A second independent solution can be obtained by 
the Bessel function of the second kind, j\j^ (f) , 

defined by 

( 154) (§) - (f)cosnT 



When n is an integer, N„(f) Is a second solution, 
but this i3 useless at f-0, for the function becomes 
Infinite at this point. Solutions can be obtained 



M 
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in terms of still other functions, the Hankel 
functions, (?) and (?) * which are merely 

linear combinations of the first and second Bessel 
functions as follows 

( 155 ) HU'lf) =t 

(156) H®(f) = 



But the Hankel functions likewise become infinite 
for f - 0 and so are excluded from use where the 
argument i3 very small. Expansion of the J and N 
functions shows a similarity to the cosine and sine 
functions and, as such, represent standing waves 
\dien used with the time factor . The H 

functions correspond to exponential functions wi th 
imaginary exponents, and, when used with the time 
factor e' AU)/t represent traveling cylindrical waves. 



The Field of Circularly Cylindrical V/ave Functions 

Every electromagnetic field within a homogeneous 
isotropic domain can be represented by -linear 
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combinations of the following functions 

(i57) j<u4 



(158) 



% 






- 6 



X.//I6 i $ 



ti 



r) 






The first equation, (157), alone, is sufficient to 
represent an electromagnetic field within finite 
domains including the axis . But at great 

distances from the source, (158) must be used, for 
as(^* r) becomes very large, the function reduces 
to a wave travelling radially outward. This can be 
more readily understood when it is noted that for 
very large arguments (f) the Hankel function is 
approximately given by 




Each elementary wave is Identified by the parameters 
n, h, k. When n is zero the field is independent of 0 
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and is symmetric about the axis, when h is zero the 
propagation is purely radial and the field becomes 
only two dimensional. The two functions, (157) and 
(158) may represent inhomogeneous plane waves as 
well as homogeneous. The planes of constant phase 
are propagated along the axis with a velocity 
y= where c* is the real part of h, but the 

amplitudes over these nlanes are functions of T 
and 0. 

The formulas for the components of the field 
vectors can be obtained from equations (147) and (l48), 
remembering that 11^ = II 0 = 0 and II Z = ^ . 

Also, the above equations we re for a generalized 
cylindrical field, while we are now restricted to 
circularly cylindrical functions and therefore, 

Ai-r t and ty'= r J (4* -Q ) I • Note also 
that, due to the exponential^ , differen- 

tiation with respect to 2 multiplies the function 
by iUA • For the transverse magnetic field, the 
results are: 




d*» r^-^Sr 

ip = (t? -jdp 
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( 161 ) 



Similarly 



(162) 



u 0) - _ M * 1 w 
n r ^cu r de 

u (0 .. M* 

”0 ~ mu ) <?r 

Hz - 0 



for the transverse electric field 



r~& _ xmu) ^ 

tr “ r de 

r~W etfr 

H e ak<m ar 

e" * 0 



-- MX 0 

ntt _ 4 . ±Pi dlP 
n e ~ r d 9 

h” = (r-A’)f 



(163) 
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When initial conditions are prescribed over given 
plane or cylindrical surfaces, a solution can be 
constructed by suoerposing these elementary wave 
functions. For fixed values of OJ (or Jfa) and h, 
the field in cylindrical coordinates may be 
represented by the equations (164) and (I 65 ) below, 
in which the 6{ v and coefficients must be determined 
from the initial conditions. 



(164) 









4 

V 






m - -«o 






co 



03 OO CO 






/>)=-«> 



(165) «•*£ 2 f 









/»=-«> 
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Propagation along a circular cylinder 

Consider a circular cylinder of infinite length, 

of radius a, of material with propagation constant 

ki Imbedded in an infinite homogeneous medium of 

propagation constant k^ . The field must be finite 

at the center and consequently the wave functions 

’within the cylinder will be constructed of Bessel 

% 

functions of the first kind. Outside the cylinder 
the field will extend to infinity and consequently 
the functions must be constructed of Hankel functions. 

In the equations below x, a ls substituted for 

and jcf is substituted for . The function 

j_ 3 abbreviated to Fn. Summations 
are in each case from n =-«> to n . A 

prime above a cylinder function indicates differentiation 
with respect to the argument (X . The super- 
scripts 1 or e over the coefficients and 

indicate points interior and exterior to the cylinder 
respectively. 






r„ 



/K 



( 165 ) &' = - 2 + 



2 1 JtMaar*. 
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(167) He - 2 ^(w) #*- ^.JnCXirfb* 



Fa 



Oi 



fC- 2 



E r e = 2 

( 168 ) E d e = ~Z [|f r H * M ^ f ^ff ^ M #] ^ 

E* = 2 lh^M<]/* 




( 169 ) h : = 2 [^ «0 k 



•/K 



H* = 2 CH "<**■)«]&> 
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The coefficients of the expansions and the 
propagation factor h must next be determined. At 
the boundary where r is equal to the constant radius 
of the cylinder, a, the tangential components of the 
field are continuous, which condition supplies a 
relation between the coefficients. From the 
continuity of the tangential components of E, the 
following relations are found, where the substitutions 
have been -made X,a -P , X a a-AT 



From the tangential components of H, similar 
relations are found as follows: 



(170) 




On) J,(r)ai ‘ H iw)af 



(172) 




H °J - ^75 $(y)hS 
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d73) = HS’WAf 



The above relations constitute a homogeneous 
system of linear equations satisfied by the four 
coefficients CtZ J b„ t Q* t . There will be 

a nontrivial solution only if the rank of the matrix 
is less than four, which will be true if its 
determinant is zero. The propagation factor is 
determined by the condition that the determinant 
of (170) through (173) shall vanish. 

Expansion of this determinant leads to the 
transcendental equation 



(174) 



six (e! _ 



J 

WR 



T& (Pj 

«,P 



-lilAll 




The roots of (174) are the allowed values of the 
propagation factor, h, and so determine the 
characteristic or natural modes of propagation. 

Thus, there has been developed an eouation 
whose solution enables calculation of natural modes 
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and the propagation constants for long rods. Equation 
(149) expresses the true result of Maxwell in an 
infinitely long rod Imbedded. in a perfectly 
homogeneous medium. Equations (157) and (158) are 
valid solutions to this equation. By use of the 
continuity of the tangential components of S and H 
at the boundary of the rod, it was possible to use 
(157) and (158) to find relations which enable one 
to solve simultaneously for the propagation constants. 

Solutions to (174) have been obtained and it 
9 10 

has been shown * that they verify experimental 
results within the accuracy of experimental procedure. 
In the next chapter, a nomogram will be presented 
which provides solutions for the case of dielectric 
rods in air or vacuum. 

It should be emphasized that the preceding 
development assumes that the rod is long enough 
to neglect boundary conditions at the ends, and that 
the radius, a, is not microscopic as far as the 
external wave is concerned. The question a3 to what 
size rod Is required to justify the assumption that 
the radius is not microscopic is most pertinent. 

If the radius is small compared to a wavelength, 

it is probably in error to assume that Hankel functions 
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are required for all waves exterior to the rod 
while Bessel functions of the first kind suffice 
for internal waves. 

Having obtained the propagation constants, 
available published literature indicates that all 
authors proceed to determine radiation characteristics 
on the assumption that the propagation is truly 
constant throughout the rod. The experimental 
evidence that this is not true in short rods of 
diameters of a wavelength or less has apparently 
been neglected, and available published theories 
on powers show serious defects when confronted with 
experimental evidence. For examnle, experimental 
curves 0 ’ show that power on axis is linear with 
exposed length of dielectric for short antennas. 

The published theories do not permit such a 
relation. 

A satisfactory solution to this problem from 
a theoretical standpoint cannot be completed in 
the present work. The experimental data presented 
may provide a stepping stone toward its ultimate 



achi evenent 
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CHAPTER 2. 



APPLICATION OF ELECTROMAGNETIC THEORY TO 
CIRCULAR DIELECTRIC RODS 



In Chapter 1 an equation was developed 
expressing in circular Coordinates the natural 
modes of electromagnetic waves in an infinitely long 
cylinder imbedded in an infinite homogeneous medium. 
The next step is to specialize this to the case of 
a dielectric cylinder imbedded in air or vacuum. 
Furthermore, particular solutions of (174) are 
desired, and to achieve this, the equations will be 
put in a form more suitable for use with available 
tables of cylinder functions. 

The equation (174) may be expressed in Modified 
Bessel functions (K n functions) by use of the 
relation 



(175) K,ft) f e ** tC’W 
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or 



0 ) 






T 



- 'Mi , . 



and 



(177) 




W© let the above expression be -g, by definition, 
to simplify notation. Also let 

( 170 ^ ^'nI.-X. — 

[ 7 ' } ? Jn(P) 

Substituting these relations in (174) 

(180) 

faf f £iq) = + rp) 1 ' 



Letting -iv = a, then 



( 178 ) 
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Sine© we are only concerned with dielectric rod 
and air or free space, we may consider 

(181) - M x - M. 

The constants k and h were introduced in our 
elementary harmonic solutions of the scalar wave 
equation in a cylinder, equations (143-145). Prom 
this, k was equal VcH^AS , and we may now 

set S-0 , since the conductivity of dielectric 

and air are negligible. Thus = CJ Vm.1?, ,and 

^a.= loTIZti . In the notation Introduced with 
equations ( 166-169), X, 1 and X? , 

Since X, = jf and X^^ it follows that 

(182) 

-A* + -^r = 

(183) 

Y/here the subscript on has no w been changed to 
zero to indicate the special case of free space. 
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Eliminating h from the above, the following 
is obtained 



Let us now set for the relative dielectric constant 



vrhlch will be one of our principal relationships 
to use in solving (174) for the propagation constants. 

To obtain the other principal equation, we 
substitute into (180) the value of h in terms of 
the other constants as given in ecuation (182), 
substitute for k^ and k^ their corresponding values 
as above, and substitute for aJ from (187). After 
simplification, using (186), the result is 



(184) = 



ass) e = £ 



and use the relation 




to put equation (184) in the form 



(187) 
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(IBS) 



(f f + ^)(f +?) = /ft* (fs ■* fX? + -f ) 



For a given n, to , and £, the last. two eouations 
can "be solved for p and q X'/hich define the mode 
propagated along the rod, and enable calculation 
of the other constants. 

From equations (182) and (I 83 ) some interesting 

deductions can be drawn concerning the longitudinal 

propagation factor, h. Me note from (I83) that if 

q is zero, that h is equal to oS{jut <T , and therefore 

the wave is traveling with the velocity and wavelength 

characteristic of free space. If p is zero, (182 ) 

shows that h is eoual , and the wave is 

traveling v;ith the wavelength characteristic of the 

dielectric. For a given dielectric, freauency, 

and diameter, equation (I 87 ) shows that the sum 

(p 2 t q 2 ) is a constant. It will later be shown 

that p is subject to less variation than q. If 

q is small, which is the case for small diameters, 

the K n (q) becomes very large compared to J^Cp). 

Since the wave functions outside the cylinder were 

( 1 ) 

constructed from H functions, which were later 
converted to functions, it is apparent that the 
energy outside the dielectric is large compared to 
that inside. Also for small q, the wavelength and 



\ 
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velocity characteristics are near that of free 
soace. Conversely, for large diameters, q is 
large relative to p, and K n (a) becomes very small 
compared to J n (p) . Thus the energy is contained 
mostly Inside the dielectric, with a relatively 
thin film of energy on the outside surfaces. The 
wavelength in this case approaches the characteristic 
of the dielectric. Qualitatively, we might expect 
in view of this, unless the dielectric is ideal 
and lossless, that losses Increase rapidly with 
increasing diameter of the rod* 

An important factor needed for the design 

of dielectric antennas is the ratio of free space 

of 

wavelength to the wavelength i propagation, or 
apparent index of refraction. Vfe have noted, 
where conductivity is zero, that: 

(a) the propagation constant, characteristic 
of the dielectric, , is equal 

(b) the propagation constant characteristic 
of free space, 

k 2 = k 0 - uYZJ, = ^ 

A 

(c) the longitudinal propagation factor for 
the combination is 
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It follows that 



(189) 



JL. 

A*. 




= 



where n Q is used as a symbol for the apparent 

index of refraction. In order to obtain an explicit 

expression for n in terms of the other constants 

a 

o 

in these equations, divide ^equation (180) by 
and note that 




j ? 4 

Solving for gp yields 



(190) 




Substituting for the numerator its qui valent given 
by (188), and extracting the square root, we have 
n in terms of € , p, and q, 

3 » 

’'Yla, ~ 



(191) 



(» 

(Zf‘ + Wfr 
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uhich is in a convenient form for obtaining n^ 

when p and q have bean determined. 

Curves of calculated apparent indices of 

refraction for various diameters and dielectric 

constants for the dominant mode (n * 1) are given 

in Figure 5. These are in reasonably close agreement 

1 ? 

with experimental values obtained by McKinney . 

It can be readily understood that this 

apparent index of refraction Is of fundamental 

importance in the design of a dielectric antenna. 

8 XX X3 

It has been shown 0 * * ^ that the optimum length 
of a cylindrical dielectric rod, to maximize 
radiation on axis, is that length which causes the 
wave to emerge from the end delayed by a half period 
behind the external wave. This length has been 
s hown‘ Q >^>^3 to be 



( 192 ) 



^0 - 



oJ (*%.'!) 



which is obtained as follows. Consider that the 
time for the wave phase to pass through -the dielectric, 
t^, must be greater than the time for the phase to 
sweep through the same distance in air, t a , by one 
half period, gy . 

Thus 



+ _ + = A _ A. - JL 

A d A ck at; o 
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Substituting r = \ and v n = - M- £ 2 . 

ana solving for i 0 * equation (192) is obtained. 

The above expression for the optimum length has 
been verified experimentally by several investi- 
gators^* 

In the discussion so far, the mode of the 
electromagnetic waves, which Is determined by the 
n in the preceding equations, has not been 
specified and the equations hold for any mode. 

It was shown in Chapter 1 that, if the wave is to 
be periodic in 0, that n must be an integer. For 

1 

a practical dielectric antenna, it is necessary 

to specify this further. The mode normally used 

for propagation in the wave guide and for exciting 

the dielectric is the so called dominant mode, 

which Is the mode of propagation having the lowest 

4 

cut off frequency. Shelkunoff has shown , by 
allowing q to approach zero and evaluating the 
remaining terms In the equation, that the n * 1 
mode has a cut off frequency which is theoretically 
zero. For n zero or greater than 1, it is also 
shown that the cut off is not zero. Therefore, 
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the dominant mode is n * 1 and this is the only 
mode considered for the remainder of this work. 

The method of solution of equations (187) 
and (1S8) must be either numerical or graphical, 
and both of these methods are cumbersome and 
tedious. In order to facilitate obtaining rapid 
solutions, at a small sacrifice of precision, 
nomograph charts A and B have been prepared. See 
Figures 6 and 7. 

Chart A, Figure 6, is a nomogram of the equation 



To obtain a solution, we seek values of p and q 
to satisfy both equations. It is assumed that the 
relative dielectric constant, 6 , the radius, a, and 
the free space wavelength, /U , are known. To use 
the accompanying charts, proceed as follows: 

Calculate the constant, p^ + , by slide rule, 

from the equation (187). Thus 



p 2 

p + q = constant 



and Chart B, Figure 7, is a nomogre.m of 



(Cf tgXf +$) - ( € /p + + V%) -0 




N omogram 

CHART A 



301 



f 



ii- 
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and locate this point on the center column of 
Chart A. A straight edge pivoted at this point 
will give the values of p and q allowed for the 
given dielectric and diameter. Assume a value for 
q and place the straight edge at this point on the 
left column, keeping it also at the pivot point 
of the center column. The corresponding value of 
p is read where the straight edge intersects the 
right column. 

Now look at Chart B, and mark off or visualize 
an € line, corresponding to the given value, running 
straight and parallel to the dash dot lines. Find 
the assumed value for q along this line, as indicated 
by the dashed q lines. At this point, thus located, 
read the corresponding value of p by eye inter- 
polation between the solid p linos. Compare the 
values of p as given by Chart B with that obtained 
from Chart A. If the value of p as given by Chart B 
exceeds that given by Chart A, then the assumed 
value of q was too large. Assume a lower value of 
q and again find the corresponding p values from 
Charts A and B as before, h’hen the two charts 
give identical p's for the same q, the solution 
is correct. 
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The above is a trial and error method, but 
it is easily performed and converses very quickly 
to a solution which i3 sufficiently accurate for 
many purposes. It may happen that a solution may 
be required that does not fall on the chart. It 
was impossible to make a chart which was readable 
and of reasonable size without considerable 
restriction on the range. However, the range 
selected for the accompanying charts is believed 
to be useable for most of the practical problems 
anticipated at this time. If the required solution 
does not fall on the charts, numerical methods should 
be used. In order to facilitate these, the values 
of f and g as functions of p and q, respectively, 
may be obtained from the plots of these functions, 
given in Figures 8 and 9. These plots may also be 
useful to refine the graphic solutions obtained 
from the Charts A and B, if additional accuracy is 
desired. 

Having solved the equations for the waves in 
a circular dielectric rod in air, we can now 
readily calculate the optimum length for a given 
dielectric constant, diameter , and frequency 
(equation 192). The next problem is to determine 










- •»! 

x : 



i» V) ih*' ^ mcU. MU kf«<* 




73 



what gain can be expected from the antenna. There 

is no accepted method for calculating the gain, 

11 15 

although several theories have been proposed * . 

It is for this reason that vie turn to experiment 
to obtain the gain to be expected for the given 
conditions. The next chapter will describe the 
equipment and procedure used to obtain the needed 
experimental data-. After describing these 
measurements, the data will be discussed and an 
attempt made to correlate the results with Marsten's 
theory"^. 

It will be noted that the preceding theory 
has dealt exclusively with the circular cross section 
rod antennas. The case of the rectangular cross 
section is much more difficult due to the fact that 
there are discontinuities on the edges, which are 
cumbersome if not impossible to account for. Further- 
more, there is no convenient coordinate system in 
which to set down the mathematical expressions. 

No published theory of electromagnetic waves In 
'dielectric rods of rectangular cross section is 
available, and consequently for the design of dielectric 
antennas of these shapes, the experimental data 
must necessarily be used. It will later be evident 
that the performance of rectangular rods is very 
de 
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definitely different from that of circular rods. 
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CHAPTER 3. 

DESCRIPTION OF EQUIPMENT, PROCEDURE, AND EXPERIMENTAL DATA 

The objective of the experimental work was to 

determine the gain on axis for dielectric antennas 

excited in the dominant mode, using optimum length 

rods in each case. These were measured for cylindrical 

rods of varying cross section diameters, and for 

rectangular rods with varying sizes in the E dimension, 

holding the H dimension constant at the width of 

the standard 3 centimeter wave guide (0.9 inches). 

A wavelength ( A e ) of 3.2 centimeters, equivalent 

o 

to a freouency of 9.37 x Kr cycles per second, was 
used for the reason that the components were of 

convenient size and the equipment for use at this 

i 

frequency was available. The eaulpment lisoBho&n in 

r 

Figures 10-14. Two dielectric materials were used 
for each antenna: Polystyrene ( £ - ) and Teflon 

( € =■ J./o ). 

A standard Klystron oscillator was used for a 
transmitter, but square wave modulation was imposed 
on the carrier frequency to facilitate detection in 
the receiver. The transmitter operated from a 
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regulated power 3upply which was fed from the 60 
cycle 115 volt power lines through a voltage 
regulating transformer. The receiver consisted of 
a crystal detector, amplifier, and voltmeter. As 
a check on receiver operation, all measurements 
were repeated using a Pickard and Burns Bolometei? 
Amplifier, model 100. Results were in substantial 
agreement with those reported and will not be 
repeated. 

The various sizes of dielectrics were fitted 
into sections of wave guide, made especially for 
the purpose, tapered down to standard size and 
fitted wi th a standard choke. See photographs. The 
part of the dielectric which fits into the guide was 
tapered in each case, in order to minimize reflections 
and reduce standing waves. The elimination or 
correction for the standing waves proved to be the 
most difficult problem to overcome. The taper 
reduced the voltage standing wave ratio (VSWR) in 
most cases to below 2.0. In every case this was 
carefully measured and a correction applied to obtain 
the actual transmitted power, using the well known formula 
(193) 

fi. _ 4r 
P ^ ~ (y h ) 2 
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where P^ Is the transmitted, oower 
is the incident power 

r is the voltage standing wave ratio (VSWR) 

Several methods were used to make the measurements. 
To begin with, the dielectric was used as the 
transmitting antenna, since receiver power was 
inadequate to accurately measure the VSWR. The VSWR 
was measured at the transmitter using a standard 
standing wave indicator. Later it was found more 
precise to use a calibrated attenuator together 
with a slotted section, crystal detector, amplifier, 
and voltmeter, attenuating as necessary to equate 
the voltage from the high voltage point to the low 
voltage point. The method by which the final 
measurements were made used this equipment, employing 
an attenuator which had been accurately calibrated 
by the National Bureau of Standards. 

As a further check on the accuracy of the 
measured gains, the various dielectric antennas 
were each put at the receiver with a special matching 
section, and the measurements repeated. For each 
reading, the matching section was adjusted to 
maximize readings by the reduction of the standing waves, 
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Results for the dielectric antennas up to one 
wavelength in diameter or in the E dimension v;ere 
in agreement with those measured by the nrevious 
two methods. Gains measured for larger antennas 
were somewhat less, which can be accounted for by 
the fact that standing waves still existed to an 
appreciable amount despite the matching section. 

Some slight difficulty was encountered with 
reflections within the room. These were reduced 
as much as possible by the strategic use of an 
absorbing screen. As a check on errors which might 
be introduced due to reflections, the distance from 
the transmitter to receiver, normally of the order 
of one hundred wavelengths, was changed in nuarter 
wavelength increments, noting power detected at 
each range. If appreciable reflections were 
interfering, there should be alternate addition 
and cancellation of power observed. It was found 
that detected power was not appreciably different 
for these positions a quarter wavelength apart, 
and it was concluded that any errors from reflections 
were of little significance. 

Before measuring the gains it was necessary 
first to find the optimum length of the dielectric. 
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This was done experimentally by observing relative 
power transmitted, first with the dielectric almost 
flush with the mouth of the guide, and then out i 
inch and so on with longer lengths using increments 
of about y X c . Near the peak of the transmission 
smaller Increments were used to obtain the exact 
optimum length. Having the relative power trans- 
mitted at optimum length, the test antenna was then 
Immediately replaced by a precalibrated standard 
horn and the powers compared. This enabled 
calculation of the relative gain of the test antenna 
with respect to the standard horn. 

The gain of the standard was carefully 

measured in order to be able to determine the 

absolute gain; that is, the ratio of the peak intensity 

of the given antenna to the intensity of an isotropic 

radiator when both are radiating the same total 

power. In terms of the effective area, f\ e , the 

oo 23 

gain, G, has been shown * to be 

U94) G ‘ 

The effective area is the ratio of the power 
received at the antenna to the power per unit area 
in the incident wave. 
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The absolute gain of the standard horn was 
measured by several of the methods which have been 
developed and published'*’ One of these^ will 
be briefly described here, which makes use of two 
identical antennas. The standard horn ’.was attached 
to the transmitter and an identical one used on the 
receiver. Under these conditions it is easily shown^ 
that 




where P R is the received power, P T is the transmitted 
power, and d is the distance between the transmitting 
and receiving antennas. To determine the ratio of 
received power to total transmitted power, the 
transmitter and receiver were coupled directly 
together, except for a calibrated attenuator between. 
Vf T lth a large but accurately measured amount of 
attenuation inserted, a particular receiver voltage 
was noted. Then the two antennas were attached 
respectively to the transmitter and receiver and 
separated by an accurately measured distance 
(several feet), and the attenuation of the transmitter 
reduced until the same receiver voltage was obtained. 
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The difference in the amount of attenuation in the 
two cases gave the ratio of P^ to P^, in decibels, 
from which G could easily be calculated. 



Presentation of Data 

Three alternatives in the presentation of 
experimental data have been considered, and these 
will be briefly mentioned to show the reasons for 
the selection of the one given. 

(a) Report the relative gain of the dielectric 
antenna compared to the gain of the feeding v/ave 
guide with the dielectric removed. This would 
permit a direct comparison of the power on axis 
with and without the dielectric which would show 
clearly the effect of the dielectric in directing 
the beam. However this method was excluded, because 
the experimental data would not be reproducible 
unless the feed structures were identical. Gain 

of open mouth feeds are not precisely predictable 
and there would be no way to obtain absolute 
gain or effective area from these data. 

(b) Report the absolute gain of the dielectric 
antenna; that is the cower on axis relative to an 
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isotropic oscillator emitting the same total power. 

This is better, but here again variations in the 
metallic feed structures would cause variation in 
the absolute gain. Furthermore, as shown by 
eauatlon (194) above, the scales would have an 
artificial multiplier of 4TT. 

(c) Report the effective cross section area 
relative to the physical cross section area. This 
is the method chosen because: (l) this presentation 
eliminates the undesired effect of variation in 
feed structures; (2) as cross section areas are 
allo\/ed to grow, the effective cross section areas 
also grow, and these tend to become eaual. In other 
words, the captured energy tends toward that 
physically intercepted by this area. Therefore, 
a plot of effective area versus physical area will 
verge tcard a straight line with unit slope passing 
through the origin regardless of the dielectric 
used; and ( 3 ) an extremely useful piece of information 
concerning any antenna is its effective cross 
section area. From this figure it is a simple and 
direct calculation to obtain the total energy 
gathered by the antenna in a field of a given strength. 

This third method of presentation, In addition 
to providing the engineer the very useful information 
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of maximum capture area for various dielectrics 
of various sizes,’ Illustrates the variation of 
capture area from the actual area in the sizes 
of the order of a wavelength, and shows how these 
should converge to the 45° line as dimensions 
increase toward optical sizes. 



Choice of Scales 

Consideration was given to reporting dimensions 
directly in centimeters or inches. This would 
enable an engineer, employing vravelengths the same 
as used' here ( 3.2 cm) to use the data directly. 

However, this was rejected because the value of 
these data is not limited to any particular wavelength. 
Dimensions are therefore given in non dimensional 
form; that is, lengths in wavelengths and cross 
section areas in square wavelengths. Although this 
requires conversion of practical units to this system 
to enable usage of the data, it has the advantage 
that the data are applicable to a very wide band 
of wavelengths. 

A plot of effective area versus actual cross 
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section area for the cylindrical dielectric antennas 
measured is given in Figure 15. Figure 16 is the 
same plot for rectangular cross section antennas. 
Figures 17 and 18 are also enclosed for completeness 
giving the absolute gains of cylinders and 
rectangular cross sections, respectively, versus 
cross section area. 

It will be noted that measurements were not 
obtained for cylinders of large enough diameter 
to show convergence of the effective area toward 
the actual area. In the larger diameter circular 
cylinders, the mode tended to shift to the /n = o 
mode which gives a null on axis. No way to prevent 
this was devised. There was also a practical limit 
to the measurements on the low side, with the time 
and facilities available. For small cross section 
rods the optimum length becomes unmanageably long, 
and the thin rods tend to sag. Furthermore, these 
thin long rods require a long range in order to 
put the field meter in the far zone. 

The present technique does not permit extension 
of the range of measurement for circular rods, 
and the possibility of existence of TM modes in 
these dielectrics is a disturbing factor. The 
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apparent irregularity of the polystyrene rods may 
be caused by this phenomenon. In order to reduce 
the probability of spurious modes in the exciting 
radiator for the larger areas, conical horns were 
also tested with polystyrene cylinders, v/ith the 
following results 





Horn 


Antenna 


Horn diameter 


Physical Area 


Capture Area 


o-3-tt 

-X 


2.5 a* 


3.78 a: 


3i" 


5.25 A? 


5.24 A? 



Here the trend of antenna effective area toward 
actual area is noted. The horns were far from 
optimum from the point of view of radiation pattern, 
however, and the lens antennas were also probably 
not optimum. Therefore, convergence of optimum 
cylinders toward the expected asymptote would not 
be closely realized until the experiment reached 
much larger areas than the above figures indicate. 

In the rectangular rods, no tendency of the 
mode to shift was observed, but other practical 
difficulties prevented extension of the range of 
sizes. In order to prevent cut off in the antenna 
feed, w hen using small sizes, the H dimension was 
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set at a constant value of .71 /\ 0 . Small aroas, 
then, would have required very long and thin sections 
which would have caused the 3ame difficulties mentioned 
above for circular rods. The larger cross section 
areas, where the E dimension was large compared 
to the constant H dimension, gave radiation patterns 
which were relatively narrow in the plane parallel 
to the S dimension, and relatively broad in the 
plane parallel to the H dimension, thus giving 
reduced gain on axis. It was found that if the 
cross section area had been made more nearly souare, 
a radiation pattern more nearly circular would have 
been produced, and thereby increasing gain on 
axis. The gain on axis, being a function of the 
radiation pattern, then, is dependent on the shape 
of the cross section. For small cross sections this 
effect is small, but for cross sections of a wave- 
length or larger the effect is very pronounced. 

A way to avoid this would be to properly shape 
the front surface of the dielectric, or to taper 
the rod. Neither of these were attempted here. 

The end surfaces of all the antennas tested we re 
plane. 

It was evident that the larger rectangular 



92 



dielectrics would give higher oowers on axis if the 
cross section was approximately square. Therefore, 
as the wavelength region is approached in Figure 16, 
the larger physical cross section antennas appear 
at a disadvantage. Perhaps a better comparison 
in this region would have been to compare the antenna 
capture area to the capture area measured for the 
same metallic feed with the dielectric flush with 
the metallic aperture (zero length antenna). This 
method of reporting would have brought the larger 
areas more in accord v r ith the following table, 
measured only for polystyrene. 



Feed Horn 


Physical Area 


Capture 

Area 


Capture Area 


Dimensions 


in Xo 


of 


of Antenna 



Feed Horn 



.9 


x 1.8” 


1.02 


0.64 


1.03 




x It" 


1.42 


1.0 


1.7 


2V 


x 2 % ” 


4.76 


4.75 


5.3 



This table shows the trend, predicted 
previously, of the capture area of the optimum 
length antenna to approach the physical cross 
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section as the latter becomes many wave lengths in 
each direction. The shaping of these larger 
antennas was only attempted approximately and 
somewhat larger capture areas would be expected 
if optimum shapes had been employed. 
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CHAPTER 4. 

DISCUSSION OF RESULTS AND ATTEMPTED CORRELATION WITH 

EXISTING THEORIES 



The curves of effective area versus cross section 
area. Figures 15 and 16 of Chapter 3, show that 
considerable increases in effective areas over 
actual areas are obtainable if the physical areas 
are small (less than a square wavelength). These 
effective areas are quite high for small physical 
areas and reduce toward the physical areas as sizes 
become larger. The trend is apparent for the 
rectangular dielectrics, but could not be demonstrated 
with circular cross section because of shifting 
of the mode within the dielectric desoite excitation 
by the ft » 1 mode. 

As previously pointed out, the design of the 
dielectric antenna above a wavelength in cross 
section area is still not understood. The gain on 
axis, on which the experimental results are based, 
is a function of the shape of the total pattern. 



95 . 



The gain is measured by sampling this pattern on 
the prolonged axis of the dielectric. This soread 
in pattern is particularly apparent in these data 
in the sizes immediately above a square wavelength. 
In this region ray optics play an important role; 
and yet customary treatments by ray optics fail 
because of the considerable phase gradients due 
to the adjacent metallic structures supporting the 
radiator. In fact, such a treatment would involve 
serious difficulties. In the present project, the 
problem would be further complicated by the absence 
of an accurate means to determine the complete v;ave 
characteristics of each of the various feed horns 
used. Proper lens shaping would necessarily have 
to be based on a knowledge of the rhase front of 
the horn to be lensed. However, to show the trend 
beyond the one wavelength region toward the exnected 
45° asymptote a few measurements were attempted 
in which the face of the dielectric was shaped by 
means of trying many combinations of various si„e 
blocks so as to obtain a stepped profile. The 
best gains were found in these larger shapes when 
the orofiles resembled ootical lenses in physical 
shape, and the larger the area the more nearly the 
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true optical chape is reauired to optimize gain. 

This is in contrast to smaller cross section areas, 
where shaping of the front- face had little effect, 
and experimental measurements were very consistent. 
Evidently the phase gradients within the smaller 
dielectrics play a much less significant role. 

Because of the fact that these data are based 
on unoptimized profiles, it is apparent that gains of 
optimized shapes would be higher than here reported, 
and convergence to the 45° asymptote would occur much 
further out. It is expected that when better theories 
are available for the design of lenses in the wavelength 
region, curves similar to those of Figure 16 could be 
expected to converge monotonically to a line of unit 
slope, and to reach the asymptote only after dimensions 
are such that the areas are many square wavelengths. 

One experiment was attempted some years ago to 

"lens" a 30 inch parabolic aperture on X band, and it 

was not found possible to detect any increase in 

effective area. As the aperture grows, optical lenses 

are obtained, where effective area is strictly proportional 

25 

to actual area. In a recent report of experiments 
sponsored by the U.S. Array Signal Corps, tests were 
revealed on a spherical Lucite lens of cross 
section area 985 square wavelengths. Gains 
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reported, after correction for the losses occurring 
in Lucite, indicate an effective area of 635 square 
wavelengths. The apparent discrepancy is due to 
the fairly narrow illumination used in this 
experiment. The illumination duplicated the usual 
radar dish illuminations and resulted in the 
expected aperture utilization factor of 65 ^ 
commonly employed in radar dish antennas. This shows 
that effective areas actually do converge to the 
nhysical areas utilized in large sizes. 

Converging lenses have been made for the 
wavelengths considered herein. While these appear 
to give a tremendous power at their designed focal 
point, a more thorough discussion would show that 
they merely operate as converging devices, obeying 
to a first approximation the half power law 




= ^ power angle of radiation 



They do not increase the effective area of the 
aperture . 



Karsten' s Theory 



A very excellent theoretical treatment of the 
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dielectric antenna of circular cross section is 



by the U.S. Naval Research Laboratory. A circular 
dielectric rod is excited in any single mode, and 
the equations used for the waves are equivalent to 
those given in Chapters 1 and 2. Prom these are 
calculated equivalent electric and magnetic current 
distributions over the surfaces, in accordance with 
Schelkunoff ' s Soui valence Theorems . Then the 
radiation field is computed using the method of 
Schelkunoff, first for the field arising from the 
equivalent sources at the sides of the rod, 
and then for the equivalent sources at the end of 
the rod, . Certain factors common to both 

expressions are dropped, since only the relative 
field is of interest, and the final expressions 

for the © components are given as follows 



given by Marsten in a recent report^ published 





( 196 ) 



r-(2) 



(196) continued 

j 2 • JX- 

= - cos^fl -~r~ j e 



77 ( 7 ^ $ cose +to.q i (D-J# } (p) + 




The expressions for the other components of the 
field will not be repeated, since only the power 
on axis i3 of interest here. 

The notation is the same as given herein, 
except as follows : 



j = X ~ 777 
$ ~ ^ 771a “ 

J? — - length of rod. Since optimum length (X) 

A 

is here considered, substitute 

X * t (x Co * d ' -£) 



On axis Cos# = 1 and for optimum length X = ~ z 

L = $ ft Stf Q On axi s L~ 0 

A and B are functions which are always finite 



and non-sero 
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S = 




which uoon simplification and substitution from 
equations (187) and (194) becomes in our notation 



<5 = 



1 df +fi L 

T 




By actually substituting numbers, it is found 
that 8 has a plus value for all dimensions under 
consideration, and so the positive sign should be 
used after extracting square root. 

Assuming n =1, & = o ) £ = jfcTf) > 
the above equations reduce to 



(197) 

£ 



(0 

e 



j /_ 

TT tfoc I 



J,(r)(V+ 1) l 



(198) 

E e = - t (jS~) p -J.fr) C * £ +4t e J /£ 
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A plot of effective area for optimum length computed 
as above Is given in Figure 20, where the abscissa 
is taken as the physical area. 

The gains calculated by Marsten's theory show 
the same general trend as the experimental data 
of Chapter 3 show. The principal exceptions are, 
first, that all are much higher in value than could 
be found experimentally, and, second, it is noted 
that the theory does not hold at all when dimensions 
are small enough so that n a approaches 1. Equation 
(202), the expression for the contribution from the 
sides, has a denominator, which approaches 

zero as n approaches 1 with decreasing diameters. 

The term in the numerator also approaches zero 

as the diameter goes to zero, but this occurs much 
more slowly. The result is that powers on axis as 
indicated become extremely high for small rods. These 
enormous powers do not conform to experimental 
observations. 

The question of identical phases for the 
radiation from the side walls and from the ends 
does not appear to be correct from a physical viewpoint. 
The phase within an optimum length rod emerges one 
half period later than the wave travelling the same 
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distance in free space. Along the rod the phase 
relative to the end must certainly lead by something 
less than a half period but more than zero. If, 
in the theory, some adjustment to these ohases 
were made, it is possible that the theory might 
fit the experimental findings more closely. 

Another interesting point is noted from 
equations (196) and (197). As the diameter becomes 
very large p approaches zero. ( P- a V.tf-A 

and at large diameters h approaches k-^) . If n =1, 
then J n (p) =0, and the right sides of both equations 

the 

vanish. V/e must conclude from this, if f\ theory 
holds for large sizes, that no radiation could occur 
from an antenna of very large diameter if the 
optimum length is determined as calculated above. 



Wilkes 1 Theory 

An explanation for the peculiar behavior of 

gains in the region of areas of a fraction of a 

square wavelength is given in a theory proposed by 
24 

Wilkes on the Interaction of Radiant Energy. It 
has been shown experimentally that the effective 
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areas in this region are considerably greater than 
the physical cross section areas, although powers 
are definitely limited. Presently oubllshed theories 
on dielectric antennas of circular cross section 
would allow an unlimited growth of capture area 
with diminishing cross section areas. 

Wilkes, hoxvever, has shown experimentally that 
the power growth is in reality not unlimited, but 
that the capture area for small cross sections is 
limited. A simplified qualitative outline of his 
approach is given below to indicate the method by 
which he arrives at a solution of this hitherto 
unaccountable phenomenon. 

As the propagation down a very thin dielectric 
rod is a plane wave, the phase and group or photon 
velocities are assumed by him to be Identical. By 
an ingenious technicue, he has been able to measure 
the phase velocities in extremely short lengths of 
dielectric and to compute the velocities down to 
a zero length. In these thin dielectrics the 
energy starts at its critical velocity, for 

very short lengths and approaches c for increasing 
lengths. Knowing this rate of change of velocity 
it is possible to conceive of the photons following 
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describable ballistic trajectories. This, however, 
requires the concept of an applied force to photons 
of finite mass which exhibit properties of Inertia. 
Since the photon mass is at present unknown, except 
for de Broglie's controversial findings^ , direct 
use of ballistic principles to calculate this 
force is not possible. Nevertheless, with this 
concept in mind, a consideration of the nature 
of the forces brought to bear on the photon reveal 
a behavior closely compatible with experimental 
observations. 

Postulating the existence of forces between 
the internal and external fields, which we have 
just noted, it follov;s that there is a transfer 
of work or energy from one of these fields to the 
other. The medium of energy exchange is photons 
of a freouency determined by the excitation. 

Wilkes has found experimental evidence of 
the existence of these forces; in fact his 
calculations and experiments support the existence 
of the force as composed of two parts. One part 
is a pressure action proportional to the difference 
in velocities of the two ‘fields. The other is a 
function of the total phase lag and appears 
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gravitational in nature. The magnitude of this 
action is found ecual to the radiant pressure and 
parallels in this respect Mandelker’s theory^' of 
gravitation. 

Prom Quantum Electrodynamics, Wilkes finds 
that the quanta in the accelerated dielectric 
fields become negative and present holes. He then 
allows these holes to be filled in to the limit 
allowed by the statistical laws of a degenerate 
gas, and by means of these calculations arrives 
at the two saturation roints shown on Figure 21 as 
the lowest cross sections reported. These points 
have been confirmed by his experiments. The other 
points of Figure 21, namely those for cross sections 
greater than 0.2 square wavelengths, are a 
repetition of the data obtained in Chapter 3. 
Obviously, if the physical area drops to zero, the 
capture area will also go to zero, and this 
seemingly trivial statement permits the completion 
of the curves of Figure 21. 
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Conclusion 

It is believed that no theory based purely on 
Maxwell's eouations will suffice to explain the 
peculiar behavior of gains in the very small cross 
section region. As stated early in this naper. 
Maxwell's equations hold with certainty only for 
large scale elec troma.gne tic phenomena, and when 
dimensions are reduced to fractions of a wavelength, 
they cannot be relied upon. The true explanation of 
very small dielectric antennas must come from quantum 
electrodynamics, as proposed by V/ilkes, and the 
data obtained in Chapter 3 appear to be in close 
agreement with this theory. 

As dimensions grow larger, effective areas 
of dielectric antennas, if shaped for optimum gain, 
remain greater than the physical cross sections, but 
will tend monotonically toward these areas. After 
areas have reached several square wavelengths, 
effective area of optimum shaped dielectrics can be 
closely approximated by the actual physical area. 
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GLOSSARY OF SYMBOLS 



E Electric Field Intensity vector 

H Magnetic Field Intensity vector 

D Electric Displacement vector 

B Magnetic Induction vector 

J Current Density vector 

j> Charge density at any point. Also used as a 
general argument of a Bessel function 

P Electric polarization vector 

M Magnetic polarization vector 

Electric Inductive capacity of free space 
defined by D= £>£ 

Magnetic Inductive capacity of free space 
defined by 6=u/ 0 H 

in Magnetic susceptibility 

Dd e Electric susceptibility P = D 

6 Relative dielectric constant - -I 6 

Co 

K*n Permeability - 

A Vector Potential of the Field 

<j) Scalar Potential of the Field 

II Hertz vectorj-Slectrlc Polarization Potential 
II* Hertz vector. Magnetic Polarization Potential 
& Conductivity 
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IV Surface density of charge, or charge per unit area 
I Conduction current 

i Square root of minus one 

h Longitudinal propagation factor “ 

t Time 

e Base of natural logarithms 

-^(f) Bessel Function of order n 

f An argument of a Bessel function 

7 m Circular cylinder function, or Bessel function, 
of order n 

^(f) Bessel function of the second kind 

H^(f) Hankel function of the first kind 

K J f) Modified Bessel function 

k Intrinsic nhase constant = VmJcJ* 

m Angular velocity = frequenc y/zir 

X, A propagation constant for dielectric in the 
case of a circular cylinder in air = 

"X^. A propagation constant for air in the case of 
a circular cylinder in air l Oiu r -$} 

p A propagation constant for dielectric ~ x., a 

v A propagation constant for air in the case of 

a circular cylinder in air = a 

a Radius of a circular cylinder of dielectric 

q A propagation constant for air = 



Apparent index of refraction for dielectric rod 



Free soace v. r ave length 

Wavelength of propagation in a dielectric antenna 

Velocity of propagation in free space 

Optimum length for dielectric to maximize 
on-axis radiation 



Abbreviation for .ImM. 

PJn(p) 



Abbreviation for ■ ($1 
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